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Abstract 



Aims. We aim to place stronger lower limits on the cosmic infrared background (CIB) brightness at 24 yum, 70 fim and 160 fim and 
measure the extragalactic number counts at these wavelengths in a homogeneous way from various surveys. 

Methods. Using Spitzer legacy data over 53.6 deg 2 of various depths, we build catalogs with the same extraction method at each 
wavelength. Completeness and photometric accuracy are estimated with Monte-Carlo simulations. Number count uncertainties are 
estimated with a counts-in-cells moment method to take galaxy clustering into account. Furthermore, we use a stacking analysis to 
estimate number counts of sources not detected at 70 yum and 160 yum. This method is validated by simulations. The integration of the 
number counts gives new CIB lower limits. 

Results. Number counts reach 35 fiJy, 3.5 mly and 40 mly at 24 fim, 70 /jm, and 160 yum, respectively. We reach deeper flux densities 
of 0.38 mly at 70, and 3.1 at 160 fim with a stacking analysis. We confirm the number count turnover at 24 pm and 70 yum, and observe 
it for the first time at 160 yum at about 20 mJy, together with a power-law behavior below 10 mJy. These mid- and far-infrared counts: 
1) are homogeneously built by combining fields of different depths and sizes, providing a legacy over about three orders of magnitude 
in flux density; 2) are the deepest to date at 70 pm and 160 yum; 3) agree with previously published results in the common measured 
flux density range; 4) globally agree with the Lagache et al. (2004) model, except at 160 fim, where the model slightly overestimates 
the counts around 20 and 200 mly. 

Conclusions. These counts are integrated to estimate new CIB firm lower limits of 2.29j{5 nW.nT 2 .sr -1 , 5.4^5-1 nW.nr 2 .sr~', and 
at 24 yum, 70 yum, and 160 yum, respectively, and extrapolated to give new estimates of the CIB due to galaxies 
.sr _1 , 6.6^jjg nW.nr 2 .sr -1 , and 14. 6* 7 2 A 9 nW.nT 2 .sr~', respectively. Products (point spread function, counts, CIB 
contributions, software) are publicly available for download at http://www.ias.u-psud.fr/irgalaxies/ 

Key words. Cosmology: observations - Cosmology: diffuse radiation - Galaxies: statistics - Galaxies: evolution - Galaxies: photom- 
etry - Infrared: galaxies 



8.9!}-| nW.nr 2 .sr 
of 2.86*5$ nWm 



1. Introduction 

The extragalactic background light (EBL) is the relic emission 
of all processes of structure formation in the Universe. About 
half of this emission, called the Cosmic Infrared Background 
(CIB) is emitted in the 8-1000 pm range, and p eaks around 
150 um. It is essen t ially d u e to the star f ormation JPugetetal 



Fixsen et al 



1998; Hauseretal. 1998; Lagache et al 



sources were responsible for a very small part of the CIB. 
Between 1995 and 1998, the ISO (infrared space o bservatory) 
perfor med deeper observations of infrared galaxies. Elbaz et all 
(120021) resolved into the source more than half of the CIB 
at 15 pm. At larger wavelengths, the sensitivity a nd angular 
resolu tion was not sufficient to resolve the CIB dDole et al.l 
2001). 



iGispertet alJ l200Ct lHauser & Dwekl 120011: iKashlinskv 
120051: lLagache et al.ll2005h . 



The CIB spectral energy distribution (SED) is an im- 
portant constrai nt for the in frared g alaxies evolution 
model s (e.g. lLagache et alJ d2004l): iFranceschini et al] 
(120091): iLe Borene et all d2009l): iPearson & KhanT l2009); 
iRowa n-Robinson (20091; IValiante et all d2009h ). ft gives the 
budget of infrared emission since the first star. The distribution 
of the flux of sources responsible for this background is also a 
critical constraint. We propose to measure the level of the CIB 
and the flux distribution of the sources at 3 wavelengths (24 pm, 
70 yum and 160 pm). 



In the 1980's, the infrared astronomical satellite (IRAS) 
and COBE/DIRBE performed the first mid-infrared (MIR) and 
far-infrared (FIR) full-sky surveys. Nevertheless, the detected 



The Spitzer space telescope dWerner et al.l 120 04). launched 
in 2003, has performed deep infrared observations on wide 
fields. The mul tiband imaging photometers for Spitzer (MIPS) 
dRieke et all2004l) mapped the sky at 24 pm, 70 u m and 160yum. 
Abou t 60% of the C IB was resolved at 24 pm ( Pa povich et al.l 
l200l and at 70 u rn dFraver et al.l l2006). Because of confusion 
(Dol e et al] l2003l) , only about 7% were resolved at 160 pm 
dDole et aljr2004l) . iDole et alJ d2006l) managed to resolve most 
of the 70 pm and 160 pm by stacking 24 pm sources. 

The cold mission of Spitzer is over, and lots of data are 
now public. We present extragalactic number counts built 
homogeneously by combining deep and wide fields. The large 
sky surface used significantly reduces uncertainties on number 
counts. In order to obtain very deep FIR number counts, we 
used a stacking analysis and estimate the level of the CIB in the 
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Field name 


Surface area 


80% completeness flux 


Scaling factor 


24 (im 


70 fim 


160 fim 


24 fim 


70 fim 


160 fim 


24 fim 


70 fim 


160 fim 


deg- 


fijy 


rr 


Jy 




FIDEL eCDFS 


U.zj 


n 1 o 




60. 


4.6 


- 


1.1)13/ 


i 
1 




FIDEL EGS 






U.JO 


76. 




45. 


1.0157 




0.93 


COSMOS 


2.73 


2.41 


2.58 


OA 

vo. 


/.y 


40. 


1 


0.92 


0.96 


o W 1 Kc Lrl 


10.04 


11.88 


11.10 


282. 


25.4 


92. 


1.0509 


1.10 


0.93 


SWTRF FN1 


9.98 


9.98 


9.30 


261. 


24.7 


94. 


1.0509 


1.10 


0.93 


SWIRE EN2 


5.36 


5.34 


4.98 


267. 


26.0 


90. 


1.0509 


1.10 


0.98 


SWIRE ESI 


7.45 


7.43 


6.71 


411. 


36.4 


130. 


1.0509 


1.10 


0.98 


SWIRE CDFS 


8.42 


8.28 


7.87 


281. 


24.7 


88. 


1.0509 


1.10 


0.98 


SWIRE XMM 


8.93 






351. 






1.0509 






Total 


53.55 


45.51 


42.91 





Table 1. Size, 80% completeness flux density and calibration scaling factor (see Sect. 12. U of the used fields. Some fields are not 
used at all wavelengths. 



three MIPS bands with them. 



2. Data, source extraction and photometry 

2.1. Data 

We took the public Spitzer mosaics^ from different observation 
programs: the GOODS/FIDEL (PI: M. Dickinson), COSMOS 
(PI: D. Sanders) and SWIRE (PI: C. Lonsdale). We used only 
the central part of each field, which was defined by a cut of 50% 
of the median coverage for SWIRE fields and 80 % for the other. 
The total area covers 53.6 deg 2 , 45.5 deg 2 , 42.9 deg 2 at 24 /im, 
70 fim and 160 fim respectively. The surface of the deep fields 
(FIDEL, COSMOS) is about 3.5 deg 2 . Some fields were not 
used at all wavelengths for different reasons: There is no public 
release of FIDEL CDFS data at 160 fim; the pixels of the EGS 
70 fim are not square; XMM is not observed at 70 and 160 fim. 
TableQ~]summarises the field names, sizes and completenesses. 

In 2006, new calibr a tion factors were a d opted for M IPS 
dEneelbracht et al.1 120071 iGordon et all 120071: IStansberrv et al.l 
120071) . The conversion factor from instrumental unit to MJy/sr is 
0.0454 (resp. 702 and 41.7) at 24 fim (resp. 70 fim and 160 fim). 
The COSMOS G03 and SWIRE (released 22 Dec. 2006) mo- 
saics were generated with the new calibration. The FIDEL mo- 
saics were obtained with other factors at 24 fim and 160 fim 
(resp. 0.0447 and 44.7). The 70 fim and 160 fim COSMOS mo- 
saics were color corrected (see Sect. 12.31 1. Consequently we ap- 
plied a scaling factor (see Table [TJ before the source extraction 
to each mosaic to work on a homogeneous sample of maps (new 
calibration and no color correction). 

2.2. Source extraction and photometry 

The goal is to build homogeneous number counts with well- 
controlled systematics and high statistics. However, the fields 
present various sizes and depths. We thus employed a single 
extraction method at a given wavelength, allowing the heteroge- 
neous datasets to combine in a coherent way. 



2.2.1. Mid-IR/far-IR differences 

The MIR (24 fim) and FIR (70 fim and 160 fim) maps have 
different properties: in the MIR, we obse rve lots of fain t 
blended sources; in the FIR, due to confusion dDole et al.l l2004). 
all these faint blended sources are only seen as background 
fluctuations. Consequently, we used different extraction and 
photometry methods for each wavelength. In the MIR, the 
priority is the deblendin g: accordingly we took the SExtractor 
(lBertin& Arnoutslll996l) and PSF fitting. In the FIR, we used 
efficient methods with strong background fluctuations: wavelet 
filtering, threshold detection and aperture photometry. 



2.2.2. Point spread function (PSF) 

The 24 fim empirical PSF of each field is generated with the 
IRAF (i mage reductio n and analysis facilitjQ) DAOPHOT 
package (Stetson 1987) on the 30 brightest sources of each 
map. It is normalized in a 12 arcsec radius aperture . Aperture 
correction (1.19) is computed with the S Tiny TirrQ(Krist 2006) 
theoretical PSF for a constant vS v spectrum. The difference of 
correction between a S v = v~ 2 and a v 2 spectrum is less than 
2%. So, the hypothesis on the input spectrum is not critical for 
the PSF normalization. 

At 70 fim and 160 fim, we built a single empirical PSF 
from the SW IRE fields. We used the Starfinder PSF extrac- 
tion routine dDiolaiti et al.l l2000h . which median-stacks the 
brightest non-saturated sources (100 mJy < S 70 < 10 Jy and 
300 mJy< S 160 < 1 Jy)- Previously, fainter neighboring sources 
were subtracted with a first estimation of the PSF. At 70 fim 
(resp. 160 fim), the normalization is done in a 35 arcsec (resp. 
80 arcsec) aperture, with a sky annulus between 75 arcsec 
and 125 arcsec (resp. 150 arcsec and 250 arcsec); the aperture 
correction was 1.21 (resp. 1.20). The theoretical signal in the 
sky annulus and the aperture correction were computed with 
the S Tiny Tim Spitzer PSF for a constant vS v spectrum. These 
parameters do not vary more than 5 % with the spectrum of 
sources. Pixels that were affected by the temporal median 
filtering artifact, which was sometimes present around bright 
sources, were masked prior to these operations. 



1 from the Spitzer Science Center website: 
http://data.spitzer.caltech.edu/popular/ 



2 http://iraf.noao.edu/ 

3 http://ssc.spitzer.caltech.edu/archanaly/contributed/stinytim/ 
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2.2.3. Source extraction and photometry 



At 24 fim, we detected sources with SExtractor. We chose a 
Gaussian filter (gauss_5.0_9x9.conv) and a background filter of 
the size of 64 x 64 pixels. The detection and analysis thresholds 
were tuned for each field. We performed PSF fitting photometry 
with the DAOPHOT allstar routine. This routine is very efficient 
for blended sources flux measurement. 

At 7 um and 160 um , we applied the a-trou wavelet 
filtering dStarck et al.l 1 199S>h on the maps to remove the large 
scale fluctuations (10 pixels) on which we performed the sourc e 
detection with a threshold algorithm dDole et all 1200 lL 12004). 
The threshold was tuned for each field. Photometry was done 
by aperture photometry on a non filtered map at the positions 
found on the wavelet filtered map. At 70 fim, we used 10 arcsec 
aperture radius and a 18 arcsec to 39 arcsec sky annulus. At 
160 fim, we used an aperture of 20 arcsec and a 40 arcsec to 75 
arcsec annulus. Aperture corrections were computed with the 
normalized empirical PSF: 3.22 at 70 fim and 3.60 at 160 fim. 
In order to estimate the uncertainty on this correction, aperture 
corrections were computed using five PSF built on five different 
SWIRE fields. The uncertainty is 1.5% at 70 fim and 4.5% at 
160 fim. 



1000.0 



100.0 



10.0 



Detected sources 
Spurious sources 




70 fim FIDEL eCDFS 



10 



S [mJy] 

70 J 1 



[h] 



Figure 1. Flux distribution of sources extracted from normal 
(solid line) and flipped (dash line) maps, at 70 fim in FIDEL 
eCDFS. The vertical dashed line represents the 80% complete- 
ness flux density. 



2.3. Color correction 

The MIPS calibration factors were calculated for a 10000 K 
blackbody (MIPS Data Handbook 20070). However, the galax- 
ies SED (spectral energy distribution) are different and the MIPS 
photometric bands are large (A/AA » 3). Thus, colo r corrections 
were n eeded. We used (like Shup e et alj d2008l) and lFraver et alj 
(2009)) a constant vS v spectrum at 24 /_m, 70 fim and 160 fim. 
Consequently, all fluxes were divided by 0.961, 0.918 and 0.959 
at 24 fim, 70 fim and 160 fim due to this color correction. 
Another possible convention is vS v oc y -1 , This convention is 
more relevant for the local sources at 160 fim, whose spectrum 
decreases quickly with wavelength. Nevertheless, the redshifted 
sources studied by stacking are seen at their peak of the cold 
dust emission, and their SED agrees better with the constant vS v 
convention. The difference of color correction between these 
two conventions is less than 2 %, and this choice is thus not 
critical. We consequently chose the co nstant vS v convent ion to 
more easily compa re our results with IShupe et ail d2008l) and 
iFraveret al.l d2009l) . 



3. Catalog properties 

3.1. Spurious sources 

Our statistical analysis may suffer from spurious sources. We 
have to estimate how many false detections are present in a map 
and what their flux distribution is. To do so, we built a catalog 
with the flipped map. To build this flipped map, we multiplied 
the values of the pixels of the original map by a factor of -1. 
Detection and photometry parameters were exactly the same as 
for normal catalogs. At 24 /im, there are few spurious sources 
(< 10%) in bins brighter than the 80% completeness limit flux 
density. At 70 fim and 160 fim, fluctuations of the background 
due to unresolved faint sources are responsible for spurious 
detections. Nevertheless, the ratio between detected source 

4 http://ssc.spitzer.caltech.edu/mips/dh/ 



numbers and spurious source numbers stayed reasonable (below 
0.2) down to the 80% completeness limit (see the example of 
FIDEL CDFS at 70 fim in Fig.0. 



3.2. Completeness 

The completeness is the probability to extract a source of a 
given flux. To estimate it, we added artificial sources (based on 
empirical PSF) on the initial map and looked for a detection in 
a 2 arcsec radius at 24 fim around the initial position (8 arcsec 
at 70 fim and 16 arcsec at 160 fim). This operation was done for 
different fluxes with a Monte-Carlo simulation. We chose the 
number of artificial sources in each realization in a way that they 
have less than 1 % probability to fall at a distance shorter than 2 
PSF FWHM (full width at half maximum). The completeness is 
plotted in Fig. [2] and the 80% completeness level is reported in 
Tabled] 



3.3. Photometric accuracy 

The photometric accuracy was checked with the same Monte- 
Carlo simulation. For different input fluxes, we built histograms 
of measured fluxes and computed the median and scatter of 
these distributions. At lower fluxes, fluxes are overestimated 
and errors are larger. These informations were used to estimate 
the Eddington bias (see next section). The photometric accu- 
racy at 70 fim in FIDEL CDFS is plotted as an example in Fig. [3] 

We also compared our catalogs with published catalogs. 
At 24 um, we com pared it with the GOODS C DFS catalog of 
ICharv etalJ d2004l) . and the COSMOS catalog of lLeFloc'h et all 
(2009). Their fluxes were multiplied by a corrective factor to be 
compatible with the vS Y = constant convention. Sources were 
considered to be the same if they are separated by less than 2 
arcsec. We computed the standard deviation of the distribution 
of the ratio between our and their catalogs. In a 80-120 fiJy bin 
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Figure 2. Completeness at 24 fim (left), 70 pm (center), and 160 //m (right) as a function of the source flux for all fields. The dashed 
line represents 80% completeness. 



3.4. Eddington bias 




■H-h- 



I I I I 



70 /urn FIDEL eCDFS 



10 

S [mJy] 

70 • 1 



100 



[hj 



Figure 3. Ratio between measured flux and input flux computed 
from Monte Carlo simulations at 70 pun in FIDEL eCDFS. Error 
bars represent 1 cr dispersion. The vertical dashed line represents 
the 80% completeness flux density. 




~s 1.0 



10 
S [mJy] 

70 L JJ 

Figure 4. Eddington bias: ratio between the number of detected 
sources and the number of input sources at 70 pm in FIDEL 
eCDFS. The vertical dashed line represents the 80% complete- 
ness flux density. 



in the CDFS, we found a dispersion of 19%. In a 150-250 pJy 
bin in COSMOS, we found a scatter of 13%. The offset is +3% 
with COSMOS catalog and -1% with GOODS catalog. At 70 
and 160 /an, we compared our catalogs with the COSMOS and 
SWIRE team ones. In all cases, the scatter is less than 15%, and 
the offset is less than 3%. At all wavelengths and for all fields, 
the offset is less than the calibration uncertainty. 



When sources become fainter, photometric errors increase. 
In addition, fainter sources are more numerous than brighter 
ones (in general dN/dS ~ S~ r ). Consequently, the number of 
sources in fain t bins are overestimate d. This is the classical 
Eddington bias (lEddingtonll913lll940h . The example of FIDEL 
CDFS at 70 pm is plotted in Fig. ffl 
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To correct for this effect at 70 fim and 160 fim, we estimated 
a correction factor for each flux bin. We generated an input 
flux catalog with a power-law distribution (r = 1.6 at 70 fim, 
r = 3 at 160 ^im). We took into account completeness and 
photometric errors (coming from Monte-Carlo simulations) to 
generate a mock catalog. We then computed the ratio between 
the number of mock sources found in a bin and the number of 
input sources. This task was done for all fields. This correction 
is more important for large r (at 160 fim). At 24 fim, thanks 
to the PSF fitting, the photometric error is more reduced and 
symmetrical. Less faint sources are thus placed in brighter flux 
bins. Because of this property and the low r (about 1.45), this 
correction can be ignored for 24 fim counts. 



account a difference of the reference SED: 10 000K versus 
constant vS v , and by another correction of 3% corresponding to 
the offset observed in Sect. 13.31 There is an excellent agreement 
with their work. 

The lLagache et all d2004tPI and iLe Borgne etal] (1200911*1 
generally agree well with the data, in particular on the faint 
end below 100 yuJy, an d on the positi o n of the peak around 
300 fiJy. However, the lLagache et al.l (l2004h model slightly 
underestimates ( about 10%) the counts above 200 yuJy. The 
ILe Borgne et al.l (|2009) model is flatter than the data, and agrees 
reasonably well above 600yuJy. 



4. Number counts 

4.1. Removing stars from the catalogs 

To compute extragalactic number counts at 24 //m, we removed 
the stars with the K - |24] < 2 col o r crite rion and identification 
procedure following IShupe et"aT1 (120081) . The K band mag- 
nitud es were taken from the 2MASS catalog ( Skrutskie et al.l 
2006). We ignored the star contribution at 70 fim and 160 fim, 
which is negligible (<1% in all us ed flux density bins) according 
to the DIRBE Faint Count model (lArendt et alJil998h . 



4.2. 24 /jm number counts 

We counted the number of extragalactic sources for each field 
and in each flux bin. We subtracted the number of spurious 
detections (performed on the flipped map). We divided by the 
completeness. As a next step, the counts of all fields were 
combined together with a mean weighted by field size. Actually, 
a weighting by the number of sources in each field overweighs 
the denser fields and biases the counts. Counts from a field 
were combined only if the lower end of the flux bin was larger 
then or equal to the 80% completeness. We thus reached 71 fi]y 
(71 //Jy to 90 yuJy bin) in the counts. However, to probe fainter 
flux densities, we used the data from the deepest field (FIDEL 
eCDFS) between a 50 and 80 % completeness, allowing us to 
reach 35 /Jy. 

Our number counts are plotted in Fig. [5] and are writ - 
ten in Table [2] We also plot data from i Papovich et alJ (|2004), 
IShupe et alj ( 120081) and [LeFloc'het ail (120091). and model pre - 
dict ions fromlLagache et al.l d2004l) and lLe Borgne et al. (2009). 
The lPapovich et al.ld2004l) fluxes are multiplied by a factor 1 .052 
to take into account the update in the calibration, the color cor- 
rection and the PSF. This correction of flux also implies a cor- 
rection on number counts, according to: 

where S , is the initial flux, S / is the corrected flux and c the 
corrective factor (S / = cSj). A correction of the flux thus 
corresponds to a shif t in the abscissa (fa ctor c) and in the 
ordinate (factor c 1 5 ). Papo vich et al.l (120041) do not subtract 
stars and thus overestimate counts above 10 mJy. We have a 
very good agreement with their work belo w 10 mJy. W e also 
have a very good agreement with IShupe et"ai] d2008l) . The 
LeFl oc'h et al.l d2009l) fluxes are multiplied by 1.05 to take into 



4.3. 70 fim number counts 

Counts in the flux density bins brighter than the 80% com- 
pleteness limit were obtained in the same way as at 24 fim 
(Fig. and Table [3J. In addition, they were corrected from the 
Eddington bias (c.f. Sect. 13.41 1. We reached about 4.9 mJy at 
80% completeness (4.9 to 6.8 bin). We used CDFS below 80% 
completeness limit to probe fainter flux density level. We cut 
these counts at 3.5 mJy. At this flux density, the spurious rate 
reached 50%. We used a stacking analysis to probe fainter flux 
density levels (c.f. section[5]). 

We can see breaks in the counts around 10 mJy and 20 mJy. 
These breaks appear between points built with a d ifferent set 
of fiel d s. Our coun t s agre e wi th earlier works of iDole et all 
d2004l) . iFraver et al.l d2006l) and iFraver et all d2009l) . However, 
these works suppose only a Poissonian uncertainty, which 
underestimates the error bars (see Sect. I4.51 l . Our data also agree 
well with these works. The Lagache et al. (2004) model agrees 
well with our data. The ILe Borgne et a l. (2009) model gives a 
reasonable fit, despite an excess of about 30% between 3 mJy 
and 10 mJy. 



4.4. 160 fim number counts 

The 160 fim number counts were obtained exactly in the same 
way as at 70 fim. We used COSMOS and EGS to probe counts 
below the 80% completeness limit. We reached 51 mJy at 
80% completeness (51 mJy to 66 mJy bin) and 40 mJy for 
the 50% spurious rate cut (figure [7] and table @J. We used a 
stacking analysis to probe fainter flux density levels (c.f. Sect.|5]l. 

Our counts agree w ith the earlier w orks of |Dole et aE [j2004l) 
and IFraver et all (120091) We find like IFraver et all d2009h that 
the lLagache et al.l (2004) model overestimates the counts by 
about 30% above 50 mJy (see the discus sion in Sect. 17.21 ). On 
the contrary, the ILe Borgne et al.l d2009l) model underpredicts 
the counts by about 20% between 50 mJy and 150 mJy. 



4.5. Uncertainties on number counts including clustering 

IShupe et all d2008l) showed that the SWIRE field-to-field vari- 
ance is significantly higher than the Poisson noise (by a factor of 

5 lLagache et all i2004l) model used a ACDM cosmology with 
Q A =0.73, n M =0.21 a nd h = 0.71 

6 ILe Borgne et alj j2009T) model used a ACDM cosmology with 
n A =0.7, n M =0.3 and h = 0.7 
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< s > 




S max 




C poisson 


& clustering 


G"clus.+calib. 




(in mJy) 




(in gal.Jy 


.sr ) 




deg 2 


0.040 


0.035 


0.044 


17.5 


1.0 


1.1 


1.3 


0.2 


0.050 


0.044 


0.056 


21.4 


1.0 


1.1 


1.4 


0.2 


0.064 


0.056 


0.071 


28.2 


1.2 


1.5 


1.8 


0.2 


0.081 


0.071 


0.090 


36.2 


1.5 


1.9 


2.4 


0.2 


0.102 


0.090 


0.114 


52.6 


1.3 


1 .9 


2.9 


0.6 


0.130 


0.114 


0.145 


64.1 


1.0 


1 .7 


3.1 


3.4 


0.164 


0.145 


0.184 


78.7 


1.1 


2.2 


3.8 


3.4 


0.208 


0.184 


0.233 


89.8 


1.3 


2.8 


4.5 


3.4 


0.264 


0.233 


0.295 


96.5 


1.5 


3.3 


5.1 


3.4 


0.335 


0.295 


0.374 


1 12.0 


0.8 


1.8 


4.8 


37.2 


0.424 


0.374 


0.474 


103.7 


0.6 


1.7 


4.5 


46.1 


0.538 


0.474 


0.601 


91.9 


0.6 


1.5 


4.0 


53.6 


0.681 


0.601 


0.762 


81.2 


0.6 


1.5 


3.6 


53.6 


0.863 


0.762 
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65.3 
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1.5 
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3.6 
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7.615 


14.496 


46.9 


2.3 


3.6 


4.1 


53.6 


21.045 


14.496 


27.595 


36.4 


3.3 


4.4 


4.6 


53.6 
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27.595 


52.531 


43.4 


5.9 


7.7 


7.9 


53.6 


76.265 


52.531 


100.000 


47.7 


9.9 


12.0 


12.2 


53.6 



Table 2. Differential number counts at 24 fim. o~ c iusterin<> is the uncertaint y taking into account clustering (see Sect. 14. 5I >. cr c i us .+ ca iib. 
takes into account both clustering and calibration (Engelb racht et al.ll2007l) . 
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4.894 


2073. 


264. 


309. 


342. 


0.2 


5.868 
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26.159 


36.578 


2557. 


69. 


124. 


218. 


38.1 


43.862 
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112. 


164. 


228. 
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157.720 
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215.440 


2354. 


121. 


198. 


257. 


45.5 


339.800 


215.440 


464.160 


2048. 


200. 


276. 


311. 


45.5 


732.080 


464.160 


1000.000 


2349. 


381. 


500. 


526. 


45.5 



Table 3. Differential number counts at 70 pirn, cr c lustering is the unce rtainty taking into account clustering (see Sect. [43). a- c i us . +C aiib. 
takes into account both clustering and calibration dGordon et al .2007]). 
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45.747 
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51.493 
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444.200 
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1000.000 
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Table 4. Differential number counts at 160 fim. CT ciusMiiis. i s me uncertain ty taking into account clustering (see Sect. 14.5b . o- c i us .+ C aiib. 
takes into account both clustering and calibration (Stansberry et al. 2007). 
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Figure 5. Differential number counts at 24 pm. Filled c ircle: points obtained with >80% completeness; filled diamond: points ob- 
tained with a 50% to 80% comple t eness; open triangle: iPapovich et al.l (120041) GTO number counts obtained with P SF fitting pho- 
tometr y; open square: \Shupe et al] ( l2008l) SWIRE number counts obtained with aperture p hotometry; open diam ond: \LsFloc'h et al.l 
(2009) COSMO S number counts obtaine d with PSF fitting photometry; continuous line: Lagac he et alj ((2004) model; dashed line 
and grey region: Le Borgn e et alJ j2009) model and 90% confidence region. Error bars take into accounts clustering (see section 
14. 5i and calibration uncertainties (Engelbracht et al. 2007). 



three in some flux bins). They estimated their uncertainties on The details are explained in the appendixlAl 
number counts with a field bootstrap method. We used a more 
formal method to deal with this problem. 



The uncertainties on the number counts are Poissonian 
only if sources are distributed uniformly. But, actually, the 
infrared galaxies are clustered. The uncertainties must thus be 
computed taking into account clustering. We first measured 
the source clustering as a function of the flu x density with 
the counts-in- c ells moments (c-in-c ) method (iPeeblesI 1 1980: 
ISzapudil 119981: iBlake & WalJ l2002h . We then computed the 
uncertainties knowing these clustering properties of the sources, 
the source density in the flux density bins and the field shapes. 



This statistical uncer tainty can be combined with the Spitzer 
calibration uncertainty (Engelbra cht et alJ 120071: iGordon et al.l 
12.0071 : IStan sberrv et al. 2007) to compute the total uncertainty 
on differential number counts. 
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Figure 6. Differential number counts at 70 /vm. Filled circle: points obtained with > 80% completeness; filled diamond: points 
obtained with less than 50% spurious s ources and less tha n 80% completeness;^Wec/ square: stacking number counts ( clear: FIDEL 
eCDF S, dark: COSMOS ); open square:\Dole et alJ(l2004h numb er counts in CDFS, Bootes and Marano; open fna«g/e: lFraver etafl 
(2006]) in G OODS andlFraver et all d2009l) in COSMOS: cross: iFraver et all d2006l) deduced f rom background fluctuations; con- 
tinuous line: lLagache et al.l d2004l) model: dashed line and grey region: iLe Borgne etai] d2 009) model and 90% confidence region. 
Error bars take into account clustering (see Sect. 14.51 ) and calibration uncertainties (Gordon et al. 2007|). 



5. Deeper FIR number counts using a stacking 
analysis 

5.1. Method 

The number counts derived in Sect.|4]show that down to the 80% 
completeness limit, the source surface density is 24100 deg~ 2 , 
1200 deg~ 2 , and 220 deg~ 2 at 24, 70, and 160 fim, respectively, 
i.e. 20 times (resp. 110 times) higher at 24 fim than at 70 fim 
(resp. 160 /mi). These differences can be explained by the 
angular resolution decreasing with increasing wavelength, thus 
increasing confusion, and the noise properties of the detectors. 
There are thus many 24 fim sources without detected FIR 
counterparts. If we want to probe deeper into the FIR number 
counts, we can take advantage of the information provided by 



the 24 fim data, namely the existence of infrared galaxies not 
necessarily detected in the FIR, and their positions. 

We used a stacking analysis dDole et al.ll2006l) to determine 
the FIR/MIR color as a function of the MIR flux. With this 
information, we can convert MIR counts into FIR counts. The 
stacking technique consists in piling up very faint far-infrared 
galaxies which are not detected individually, but are detected at 
24 fim. For this purpose, it makes use of the 24 fim data prior 
to tracking their undetected counterpart at 70 fim and 160 fim, 
where most of the bolometri c luminosity arises. This method 
was used bv lDole et al.l (12006). who managed to resolve the FIR 
CIB using 24 ^m s ource s positions, as wel l many other author s 
(e.g lSerieant etal] d2004l) : lDve etaTlll2006l) : IWang et al.l d2006): 
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Figure 7. Differential number counts at 160 fim. Filled circle: points obtained with > 80% completeness; filled diamond: points 
obtained with less than 50% spurious sources and less than 80% com pleteness; fi l led sq uare: stacking number counts (clear: 
FIDEL/GTO C DFS, middle : COSM OS, dark: SWIRE EN1); open s auare:\T>o\e et alJ(f2 004) number counts in CDFS and Marano; 
open triangle: | Fraver et all d2009h in COSMOS; continuous line: lLag ache et al. (2004) model; dashed line and grey region: 
iLe Borgne et al J d2009h model and 9 0% confidence region. Error bars take into account clustering (see Sect. 14.51 ) and calibration 
uncertainties IStansberr vet alJl2007h . 



Devli n etalJ (12009 ): iDve et alJ d2009); Marsden etall (|2009); stacking. 



Pasca le et alj §009)) 



To derive the 70 yum or 160 //m versus 24 pirn color, we 
stacked the FIR maps (cleaned of bright sources) at the positions 
of the 24 fj.m sources sorted by flux, and performed aperture pho- 
tometry (same parameters as in Sect.|2](. We thus get 



S FIR - f(S 24), 



(2) 



We checked that / is a smooth monotonic function, in agree- 
ment with the expectation that the color varies smoothly with 
the redshift and the galaxy emission properties. Assuming that 
the individual sources follow this relation exactly, the FIR num- 
ber counts could be deduced from 



where S fir is the average flux density in the FIR of the 
population selected at 24 yum, S 24 the average flux density at 24 
microns, and f is the function linking both quantities. We derive 
f empirically using the Sfir versus S24 relation obtained from 



dN 



dS 1 



dN 



?=/(S 2 4) 



dS 



24 



idS 1 



dS 



24 



(3) 
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In practice, the two first terms are discrete. In addition, the last 
term is computed numerically in the same S 24 bin, using the two 
neighboring flux density bins (k-1 and k+1). We finally get 



dN 



dS 



(S 



FIRk 



) = 



FIR 



dN 



dS 



24 Ik 



> FIR,k+ 1 



1 FIR,k- 1 



> 24,i+l 



-s 



(4) 



24,*- 1 



where < S fir > is measured by stacking. In reality, sources do 
not follow Eq. [2] exactly, but exhibit a scatter around this mean 
relation 



S fir = f(S 24) + cr- 



(5) 



Our method is still valid under the condition j <K 1, and we 
verify its validity with simulations (see next section). 

To obtain a better signal to noise ratio, we cleaned the 
resolved bright sources from the FIR maps prior to stacking. 
We used 8 S 24 bins per decade. We stacked a source only if 
the coverage was more than half of the median coverage of the 
map. Uncertainties on the FIR mean flux were estimated with 
a bootstrap method. Furthermore, knowing the uncertainties on 
the 24 jjm number counts and the mean S 24 fluxes, we deduced 
the uncertainties on the FIR number counts according to Eq. [4] 

At 70 /im, we used the FIDEL eCDFS (cleaned at 
S 70 > 10 mJy) and the COSMOS (cleaned at S 70 > 50 mJy) 
fields. At 160 yum, we used 160 /um the GTO CDFS (cleaned at 
S i6o > 60 mJy), the COSMOS (cleaned at S l6Q > 100 mJy) and 
the SWIRE EN1 (no clean to probe the S 160 > 20 mJy sources) 
fields. 



5.2. Validation on simulations 



We used the iFernandez-Conde et al.l d2008l) simulations^ to 
validate our metho d. These simulations are based on the 
iLagache et al.l d2004h model, and include galaxy clustering. We 
employed 20 simulated mock catalogs of a 2.9 deg 2 field each, 
containing about 870 000 mock sources each. The simulated 
maps have the same pixel size as the actual Spitzer mosaics 
(1.2", 4", and 8" at 24 //m, 70 yum, and 160 fim, resp.) and 
are convolved with our empirical PSF A constant standard 
deviation Gaussian noise was added. We applied the same 
method as for the real data to produce stacking number counts. 
At 160 /um, for bins below 15 mJy, we cleaned the sources 
brighter than 50 mJy. Figure [8] shows 160 /mi number counts 
from the mock catalogs down to 5160 = 1 mJy (diamond) and 
the number counts deduced from the stacking analysis described 
in the previous section (triangle). The error bars on the figure 
are the standard deviations of the 20 realization. There is a good 
agreement between the stacking counts and the classical counts 
(better than 15%). Nevertheless, we observed a systematic bias, 
intrinsic to the method, of about 10% in some flux density 
bins. We thus combined this 10% error with the statistical 
uncertainties to compute our error bars. We also validated the 
estimation of the statistical uncertainty in the stacking counts: 
we check that the dispersion of the counts obtained by stacking, 
coming from different realizations, was compatible with our 
estimation of statistical uncertainties. The results are the same 
at 70 yum. 
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Figure 8. Simulated number counts at 160 /-im, computed from 
stacking counts ( triangle ) and from the input mock catalog ( solid 
line). The good agreement (better than 15%) validates the stack- 
ing counts method (see Sect. 15.21 ). Twenty realizations of the 2.9 
deg 2 field maps with about 870 000 mock sources each were 
used. 



5.3. Results 

At 70 //m, the stacking number counts reach 0.38 mJy (see 
Fig. [6] and Table [5}. The last stacking point is compatible with 
the iFraver et al.l (|2006) P(D) c onstraint. The s tack ing points 
also agree very wel l with the lLagache et al.l d2004l) model. 
The iLe Bo rgne et al. (2009) model predicts slightly too many 
sources in 0.3 to 3 mJy range. The turnover around 3 mJy 
and the power law behavior of th e faint counts (S 70 < 2m Jy), 
observed by IFraver et al.l d2006l) are confirmed with a better 
accuracy. 

At 160 /mi, the stacking counts reach 3.1 mJy (see Fig. [7] 
and Table |6j. We observed for the first time a turnover at about 
20 mJy, and a power-law decrease at smaller flux densities. The 
stacking counts are lower than the lLagache et all d2004l) model 
around 20 mJy (about 30% ). Below 15 mJy, the s tacking counts 
agree with this model. The lLe Borgne et al] d2009l) model agree 
quite well with our points below 20 mJy. The results at 160 /urn 
will be discussed in Sect. 17.21 



Publicly available at http://www.ias.u-psud.fr/irgalaxies/ 



6. New lower limits and estimates of the CIB at 
24 /im, 70 //m and 160 /im 

6.1. 24 fim CIB: lower limit and estimate 

By integrating the measured 24 /um number counts between 
35 //Jy and 0.1 Jy, we can estimate a lower value of the CIB 
at this wavelength. The counts were integrated with a trapeze 
method. We estimated the uncertainty on the integral by adding 
(on 10000 realisations) a random Gaussian error to each data 
point with the cr given by the count uncertainties taking into 
account clusteri ng. We then added the 4 % calibration error of 
the instrument dEngelbracht et al.l l2007h . We found 2.26+° ® 
nW.m _2 .sr _1 . The very bright source counts (S24 > 0.1 Jy) are 
supposed to be euclidian (dN/dS = C euc iS 2 ' 5 ). We used the 
three brightest points to estimate C euc i. We found a contribution 
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Table 5. Stacking extragalactic number counts at 70 /mi. cr c i us . 
is the uncertainty taking into account clustering (see Sect. 14.51 ). 
o~cius.+caiib. takes int o account both clustering and calibration 
(Gordon et al]|2007|). 
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2944. 
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Table 6. Stacking extragalactic number counts at 160 /mi. cr c/i(s 
is the uncertainty taking into account clustering (see Sect. 14.51 ). 
Q'rius.+rniih. takes into account both clustering and calibration 



(Stansberry et al. 2007). 
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Figure 9. Cumulative contribution to the surface brightness of 
the 24 /mi CIB as a function of S 24^. The colored area rep- 
resents the 68% and 95% confidence level. The shaded area 
represents the 5 24 < 35/iJy power-law extrapolation zone (see 
Sect. I6.lt . The 4% calibration uncertainty is not represented. 
The table corresponding to this figure is available online at 
http://www.ias.u-psud.fr/irgalaxies/ 
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Table 7. Summary of CIB results found in this article. 



We might have wanted to estimate the CIB value at 24 
jum. To do so, we needed to extrapolate the number counts 
on the faint end. Below 100 /iJy, the number counts exhibit a 
power-law behavior (Fig. [5}- We assumed that this behavior (of 
the form dN/dS = Cf a miS r ) still holds below 35 /iJy. r and 
C 'faint are determined using the four faintest b ins. We found 
r - 1 45 ±0.10 (compatible with 1.5 ± 0.1 of iPapovich et al.l 
(2004)). Our new estimate of the CIB at 24 /mi due to infrared 
galaxies is thus 2.86^j |g nW.m _2 .sr _1 . The results are plotted in 
Fig. [9] We conclude that resolved sources down to S 24 = 35//Jy 
account for 80% of the CIB at this wavelength. 



6.2. 70 /im and 160 //m CIB: lower limit and estimate 

At 70 /mi and 160 /mi, the integration of the number counts 
was done in the same way as at 24 /mi, except for the stacking 
counts, which are correlated. To compute the uncertainties 
on the integral, we added (on 10000 realizations) a Gaussian 
error simultaneously to the three quantities and completely 
recomputed the associated stacking counts: 1- the mean density 
flux given by the stacking; 2- the 24 /mi number counts; 
3- the mean 24 /mi flux density. At 70 //m, and 160 /mi, 
the ca libration uncertainty is 7 % dGordon et al.1 120071) and 
12% ( Stansberry e t al.l 120071) . respectively. We estimated 
the CIB surface brightness contribution of resolved sources 
(S 70 > 3.5 mJy and S m > 40 mJy) of 3.1^ nW.m^.sr 1 and 
l.O+o } nW.nr 2 .si--'. The contribution of 5 70 > 0.38 mJy and 
5160 > 3.1 mJy is 5.4+^ nW.m^.sr-' and 8.9!|;J nW.m^.sr 1 , 
respectively. 

Below 2 mJy at 70 fim, and 10 mJy at 160 //m, the stacking 
counts are compatible with a power-law. Like at 24 /mi, we 
assumed that this behavior can be extrapolated and determined 
the law with the five faintest bins at 70 /im, and the four faintest 
at 160 /<m. We found a slope r = 1.50 ± 0.14 at 70 /mi, and 
1.61 + 0.21 at 160 /im. The slope of t he nu mber counts at 70 /im 
is compatible with the iFraver et alj d2006l) value (1.63 ± 0.34). 
The slope at 160 /im is measured for the first time. Our new 
estimate of the CIB at 70 /im and 160 /im due to infrared 
galaxies is thus 6.6+gI nW.irr 2 .sr~', and I4.6+29 n W.m _2 .sr _I , 
respectively. We conclude that resolved and stacking-studied 
populations account for 82% and 62% of the CIB at 70 /im and 
160 //m, respectively. These results are summarized in Table [7] 
andFig.[T0]and[n] 



to the CIB of 0.032+°°^ 



nW.m 2 .sr . 



Consequently, very 

bright sources extrapolation is not critical for the CIB estimation 
(1% of CIB). The contribution of S 2 4 > 35 /Jy is thus 2.29+Jj Jg 
nW.m^.sr 1 (c.f. Table|7]i. 



7. Discussion 

7.1. New lower limits of the CIB 



The estimations of CIB based on number counts ignore a 
potential diffuse infrare d emission like dust in galaxy clusters 
ilMontier & Giardll2005l) . The extrapolation of the faint source 
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Figure 10. Cumulative contribution to the surface brightness 
of the 70 fim CIB as a function of Sjoum- The colored area 
represents the 68% and 95% confidence level. The shaded ar- 
eas represent the 0.38 < 5 70 < 3.3 mJy stacking counts zone 
and the S 70 < 0.38 mJy power-law extrapolation zone (see 
Sect. 16.2b . The 7% calibration uncertainty is not represented. 
The table corresponding to this figure is available online at 
http://www.ias.u-psud.fr/irgalaxies/ 
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Figure 11. Cumulative contribution to the surface brightness 
of the 160 fim CIB as a function of Si6o,um- The colored area 
represents the 68% and 95% confidence level. The shaded ar- 
eas represent the 3.1 < S160 < 45 mJy stacking counts zone 
and the S 160 < 3.1 mJy power-law extrapolation zone (see 
Sect. 16.21 ). The 12% calibration uncertainty is not represented. 
The table corresponding to this figure is available online at 
http://www.ias.u-psud.fr/irgalaxies/ 

counts supposes no low luminosity population, like population 
III stars or faint unseen galaxies. Accordingly, this type of 
measurement can provide in principle only a lower limit. 

At 24 fim, iPapovich et all d2004 found 2.7^ 1 nW.m^.sr 1 
using the counts and the extrapolation of the faint source counts. 



We agree with this wo rk and significantl y reduced the uncertain- 
ties on this estimation. Do le et al ] (l2006l) found a contribution of 
1.93 + 0.23 nW.m^.sr -1 for S 24 > 60//Jy sources (after dividing 
their results by 1.12 to correct an aperture error in their photom- 
etry at 24 fim). Our analysis gives 2.10 ± 0.08 nW.m'lsr 1 for a 
cut at 60 //Jy, which agrees very well. Rodighiero et al. (2006) 
gave a total value of 2.6 nW.m _2 .sr _1 , without any error bar. 
ICharv et all d2004l) found 2.0 + 0.2 nW.m^.sr -1 , by integrating 
sources between 20 and 1000 //Jy (we find 2.02 ± 0.10 for the 
same interval). 

At 70 fim, using the number counts in t he ultr a deep 
GOODS-N and a P(D) analysis, iFraver et all d2006l) found 
a 570 > 0.3 mJy source contribution to the 70 fim CIB of 
5.5 ±1.1 nW.irr 2 .sr~'. Using the stacking counts, we found 
5.5 + 0.4 nW.irr 2 .sr~' for the same cut, in excellent a greement 
and with improved uncertainties. In lDole et alJ {2006), the con- 
tribution at 70 fim of the S 24 < 60 //Jy sources was computed 
with an extrapolation of the 24 fim number counts and the 70/24 
color. They found 7.1 ± 1.0 nW.m^.sr -1 , but the uncertainty on 
the extrapolation took only into account the uncertainty on the 
70/24 color and not the uncertainty on the extrapolated 24 fim 
contribution, and was thus slightly underestimated. This is in 
agrees with our estimation. 

At 160 fim they found with the same method, 
17.4 + 2.1 nW.m _2 .sr _1 (a corrective factor of 1.3 was ap- 
plied due to an error on the map pixel size). This estimation is 
a little bit higher than our estimation, and can be explained by a 
small cont ribution (of the order of 15%) of the source clustering 
dBavouzetll2008l) . 

Our results can also be compared with direct measurements 
made by absolute photometers. These methods are biased 
by the for eground mod e ling, b ut do not ignore the extended 
emission. Fixs en et a D {[998) found a CIB brightness of 
13.7+3.0 at 160 fim, in excellent agreement with our est i mation 
(14.6^' nW.n r 2 ). From the discus sion in iDole et ail d2006l) 
(Sect. 4.1), the lLagache eTaTI d2000l) DIRBE WHAM (FIRAS 
calibration) estimation at 140 fim and 240 fim of 12 nW.m^.sr -1 
and 12.2 nW.m _2 .sr _1 can be also compared with our v alue at 
160 fim. A more recent work of Odega rd et alJ d2007l) found 
25.0+6.9 and 13.6+2.5 nW.m^.sf 1 at 140 fim and 240 //m re- 
spectively (resp. 15+5.9 nW.m _2 .sr _1 and 12.7+1.6 nW.m^.sr" 1 
with t he FIRAS scale). Using ISOPHOT data, Uuvela et all 
(2009) give an estimation of the CIB surface brightness between 
150jum and 180 fim of 20.25+6.0+5.6 nW.m^.sr 1 . 

The total brightness due to infrared galaxies at 160 fim cor- 
responds to the total CIB level at this wavelength. We thus have 
probably resolved the CIB at this wavelength. Nevertheless, 
the uncertainties are relatively large, and other minor CIB 
contributors cannot be excluded. 

In addition, upper limits can be deduced indirectly from 
blazar high energy spectrum. IStecker & de Jageri d 19971) gave 
an upper lim i t of 4 nW.m~ 2 .sr _1 at 20 fim using Mkn 421. 
iRenault et alJ (l200lh found an upper limit of 4.7 nW.m _2 .sr _1 
between 5 and 15 fim with Mkn 501. This is consistent with our 
lower limit at 24 fim. 

An update of the synthetic EBL SED of IDole et all d2006h 
with the new BLAST (balloon-born e large-aperture su bmil- 
limeter telescope) lower limits from iDevlin et alJ d2009h and 
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our values is plotted in Fig. [12] The BLAST lower limits are 
obtained by stacking of the Spitzer 24 pm sources at 250 urn, 
350 pm and 500 pm dDevlin et alj|2009t iMarsden et alll2009l) . 



7.2. 160 pm number counts 

At mo st, we observed a 30% overestimation of the lLagache et al.l 
(2004) model compared to the 160 pm number counts (Sect. l4~4l 
and 15. 3 1 and Fig. O, despite good fits at other wavelengths. This 
model uses mean SEDs of galaxies sorted into two populations 
(starburst and cold), whose luminosity functions evolve sepa- 
rately with the redshift. A possible explanation of the model 
excess is a slightly too high density of local cold galaxies. By 
decreasing the density of this local population a little, the model 
might be able to better fit the 160 pm number counts without 
significantly affecting other wavelengths, especially at 70 pm 
(more sensitive to warm dust rather than cold dust), and in the 
submillimetre range (more sensitive to redshifted cold dust at 
faint flux densities for wavelengths larger than 500 pm). 

The iLe Borgne et al.l d2009i) model slightly overpredicts 
faint 160 /mi sources, probably because of the presence of too 
many galaxies at high redshift; this trend is also seen at 70 
/mi. With our number counts as new constraints, their inversion 
should give more accurate parameters. 

Herschel was successfully launched on May 14th, 2009 
(together with Planck). It will observe infrared galaxies between 
70 /mi and 500 /mi with an improved sensitivity. It will be possi- 
ble to directly observe the cold dust spectrum of high-z ULIRG 
(ultra luminous infrared galaxy) and medium-z LIRG (luminous 
infrared galaxy). PACS (Photodetectors Array Camera and 
Spectrometer) will make photometric surveys in three bands 
centred on 70 /mi, 100 /mi and 160 /mi. Herschel will allow 
us to resolve a significant fraction o f the background a t thes e 
wavelengths (Lag ache et al.l 120031: ILe B orgne et al. 2009). 
SPIRE (spectral and photometric imaging receiver) will observe 
around 250 /mi, 350 pm and 500 pm, and will be quickly 
confusion limited. In both cases, the stacking analysis will allow 
us to probe fainter flux density levels, as it is complementary to 
Spitzer and BLAST. 



8. Conclusion 

With a large sample of public Spitzer extragalactic maps, we 
built new deep, homogeneous, high-statistics number counts in 
three MIPS bands at 24 /mi, 70 /mi and 160 /mi. 

At 24 pm, the results agree with previous works. These 
counts are derived from the widest surface ever used at this 
wavelength (53.6 deg 2 ). Using these counts, we give an accurate 
estimation of the galaxy contribution to the CIB at this wave- 
length (2.86±°;}j? nW.m^.sr 1 ). 

At 70 pm, we used the stacking method to determine the 
counts below the detection limit of individual sources, by 
reaching 0.38 mJy, allowing us to probe the faint flux density 
slope of differential number counts. With this information, we 
deduced the total contribution of galaxies to the CIB at this 
wavelength {6.6^PjiW.m~ 2 .sr~ l ). 



At 160 pm, our counts reached 3 mJy with a stacking 
analysis. We exhibited for the first time the maximum in 
differential number counts around 20 mJy and the power-law 
behavior below 10 mJy. We deduced the total contribution of 
galaxies to the CIB at this wavelength (14.6+2 <j,nW.m~ 2 •■ sr ')■ 
Herschel will likely probe flux den sities down to about 10 mJy 
at this wavelength (confusion limit, Le Bor gne et al.l (l2009l) l. 

The uncertainties on the number counts used in this work 
take carefully into account the galaxy clustering, which is 
measured with the "counts-in-cells" method. 

We presented a method to build very deep number counts 
with the information provided by shorter wavelength data 
(MIPS 24 pm) and a stacking analysis. This tool could be used 
on Herschel SPIRE data with a PACS prior to probe fainter flux 
densities in the submillimetre range. 

We publicly release on the website 
http://www.ias.u-psud.fr/irgal/, the following products: PSF, 
number counts and CIB contributions. We also release a stack- 
ing library software written in IDL. 
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Figure 12. Current measurement of the extragalactic background light spectral energy distribution from 100 nm to 1mm, with the 
cosmic optical background (COB, A < 8 yum) and cosmic infrared background (CIB, A > 8 fim). Our new points at 24 /mi, 70 fim 
and 160 fim are plotted (triangle). Lower (red) triangles correspond to the CIB resolved with the number counts and stacking counts. 
Upper (blue) triangles correspond to the total extrapo lated CIB due to infrared galaxies. BLAST lower limits at 250 yu m, 350 yum 
and 50 fjm (|pevlin et al. 200^ iMarsden et alJl2009h are represented in black arrows. The FIRAS measurements of iFixsen et al.l 
(1998) between 125 yimi and 2000 fim are plotted wi th a g rey solid line, and the 1-cr confidence region with a grey dashed line. Other 
points come from different authors (see lDole et al.l d2006l) for complete details). Old MIPS points are not plotted for clarity. 
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Appendix A: Uncertainties on number counts 
including clustering 

A. 1 . Counts-in-cells moments 

We consider a clustered population with a surface density p. The 
expected number of objects in a field of the size Q is N = pQ. 
In the Poissonian case, the standard deviation around this value 
is yfN. F or a clustered distribu tion, the standard deviation o~n is 
given by dWall & Jenkinsll2003l) 

o- N = ^jy.N 2 + N. (A.l) 
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Figure A.l. Amplitude of the auto correlation as a function of 
the flux density of the sources at 24 /mi, and best power-law fit. 



y ' 



The expected value of y is given by (Peebles 1980) 



Q 2 



(A.2) 



where w(8) is the angular two points auto correlation function of 
the sources. 



A.2. Measuring source clustering as a function of flux density 

We assume the classical power law description w(6) = 
A{S,A)6 l ~ y with an index y = 1.8. So, y depends only on A 
and on the shape of the field: 



y = A(S,A) 



J field J field 1 



Q. 2 



The uncertainty on y is given by (ISzapudi| [l998): 



'NcellN 2 



(A3) 



(A.4) 



To measure A(S,A), we cut our fields in 30' x 30' square 
boxes, in which we count the number of sources and compute the 
variance in five, three and three flux density bins at 24 yum, 70 
ytzm and 160 fim. We calculate the associate A(S,A) combining 
Eq.|AT]and|A3] 



A(S,A) = x 



Q 2 



(A.5) 



The fit of A(S 24,24 fim) versus S 24 (see Fig. IA. lb gives (x 2 = 
2.67 for five points and two fitted parameters) 



A{S,2A^m) = (2.86 ± 0.29). 10" 



— 1 

ImJy / 



0.90±0.15 



(A.6) 



The measured exponent in IA.6I of 0.90 + 0.15 corresponds to 
y/2, which is the expected value in the case of a flux-limited 
survey in an Euclidean universe filled with single luminos- 
ity sources. We fix this exponent to fit A(S 70, 70/Um) and 
A(Si6o, 160jum). We find A(lm/y, 70//m) = (0.25 ± 0.08). 10" 3 
and A(lmJy, 160/im) = (0.3 + 0.03).10~ 3 . 



Figure A.2. Relative error on the number count as a function of 
the field size. We have chosen A = 0.019 and p = 38.5 deg 2 
(values for a 80-120 mJy flux density bin at 160 jum). The field 
is a square. 



A.3. Compute uncertainties due to clustering 

With this model of A(S,A) and the field shape, we compute y 
(Eq. IA.3I >. Assuming = (N to be the number of detected 
sources in a given field and flux density bin), we deduce cr(N) 
from Eq. IA.1I and consequently the error bar on the number 
counts for a single field. 

To compute the final uncertainty on the combined counts, we 
use the following relation 



°~combM ~ 



where cr„ 



(A.7) 



v is the uncertainty on the combined number 
counts, Q, the solid angle of the i-th field, and <x, <w the uncer- 

dS 

tainty on the number counts in the i-th field (given by Var(N), 
Eq.ED. 



A.4. Discussion about clustering and number count 
uncertainties 

For a clustered distribution of sources, the uncertainties on the 
number counts are driven by two quadratically combined terms 
(Eq. IA. It : a Poissonian term V^V and a clusteri ng term -^y.Af 
(see F ig. lA2b . We have N cc Q and y oc Qfr-W 2 dBlake & Walll 
2002). When the uncertainty is dominated by the Poissonian 
term (small field), the relative uncertainty is thus proportional 
to y[Qr x l 2 . When the uncertainty is dominated by the clustering 
term (large field), the relative error is proportional to Q (1 ~ y) / 4 
(fi-°- 2 for7= 1.8). 

Consequently uncertainties decrease very slowly in the the 
clustering regime. Averaging many small independent fields 
gives more accurate counts than a big field covering the same 
surface. For example, in the clustering regime, if a field of 
10 deg 2 has a relative uncertainty of 0.2, the relative uncertainty 
is 0.2/ VlO = 0.063 for the mean of ten fields of this size, and 
0.2 x 10 ° 2 = 0.126 for a single field of 100 deg 2 . Consequently, 
if one studies the counts only, many small fields give better 
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results than one very large field. But, this is not optimal if one 
studies the spatial properties of the galaxies, which requires 
large fields. 



